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Gravity theories beyond General Relativity typically predict dipolar gravitational emission by 
compact-star binaries. This emission is sourced by “sensitivity” parameters depending on the stellar 
compactness. We introduce a general formalism to calculate these parameters, and show that 
in shift-symmetric Horndeski theories stellar sensitivities and dipolar radiation vanish, provided 
that the binary’s dynamics is perturbative (i.e. the post-Newtonian formalism is applicable) and 
cosmological-expansion effects can be neglected. This allows reproducing the binary-pulsar observed 
orbital decay. 
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General Relativity (GR) is very successful at inter¬ 
preting gravity on a huge range of scales, field strengths 
and velocities. Nevertheless, evidence for Dark Matter 
and Dark Energy may be interpreted as a breakdown 
of GR on cosmological scales. Also, GR is intrinsically 
incompatible with quantum field theory, and should be 
replaced, at high energies, by a (still unknown) quan¬ 
tum theory of gravity. Given this situation, guidance 
may come from experiments, namely those testing grav¬ 
ity in regimes involving strong fields and/or relativistic 
speeds. These experiments include measurements of the 
gravitational-wave (GW) driven orbital decay of binary 
pulsars [1], and upcoming GW interferometers [2-4]. It 
is therefore crucial to analyze gravitational emission in 
theories alternative to GR. 

Modified gravity theories typically generalize GR by 
introducing extra gravitational fields non-minimally cou¬ 
pled to the metric. An example is Fierz-Jordan-Brans- 
Dicke (FJBD) gravity [5-7], see e.g. Ref. [8] for a recent 
review of more theories. Often, even theories where no 
extra helds are explicitly added [e.g. f{R) gravity] can be 
recast as GR plus extra fields by a suitable change of vari¬ 
ables. The extra fields generally introduce “fifth forces”, 
and the motion of a body free-falling in a gravitational 
field will generally depend on the body’s nature. This 
effect can be suppressed for weakly gravitating bodies by 
assuming that the extra gravitational fields do not couple 
to matter directly, i.e. the equivalence principle (EP) can 
be restored for weakly gravitating bodies (“weak EP”). 
However, for bodies with strong self-gravity, the extra 
fields will still effectively couple to matter (because they 
are non-minimally coupled to the metric, which in turn 
is coupled to matter via gravity). This coupling will be 
increasingly important as the body’s gravitational bind¬ 
ing energy - which measures the “strength” of the body’s 
self-gravity - increases. Indeed, the extra gravitational 
fields will generally affect the body’s binding energy, and 
since in relativistic theories all forms of energy gravitate, 


the body’s gravitational mass will depend on the extra 
fields via the binding energy. As such, the inertial mass 
may differ from the gravitational mass if the binding en¬ 
ergy’s contribution to the latter is important. Indeed, 
possible deviations from the “strong” EP (i.e. the uni¬ 
versality of free-fall for strongly gravitating bodies) in 
modified gravity theories are typically parametrized by 
the “sensitivities” [9] 
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i.e. the derivatives of the gravitational mass M relative 
to the theory’s extra fields Qa; while keeping the body’s 
total baryon number N and entropy E fixed. (The sen¬ 
sitivities thus measure the body’s response to changes in 
the local value of Qa-) For weakly gravitating bodies to 
obey the weak EP, it must be ~ 0. Indeed, if the 
extra fields Qa do not couple to matter directly, they 
only enter M via the binding energy, whose contribution 
is negligible if the body’s self-gravity is weak. 

The sensitivities enter both the conservative dynam¬ 
ics of binary systems (e.g. the periastron precession) and 
the dissipative one (i.e. GW emission). Their leading- 
order dissipative effect is the emission of dipolar grav¬ 
itational radiation, i.e. binaries will produce GWs (j) ^ 
G/c^ X 0[{s^^'’ — where are the bodies’ 

sensitivities. This is a —1 post-Newtonian (PN) effect, 
i.e. it is enhanced by {v/c)~'^ compared to GR {v being 
the binary’s relative velocity), if ~ 1. There¬ 

fore, knowledge of the sensitivities is crucial to verify 
the agreement between a theory and GW observations 
(binary-pulsar data or direct detections). For example, 
binary pulsars have already placed strong constraints 
on Lorentz-symmetry violations in gravity [10, 11] and 
on certain scalar-tensor theories [8, 12, 13], and even 
stronger bounds will be possible with direct GW detec¬ 
tions [14-16]. 
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Here, we generalize previous work in scalar-tensor 
theories [17, 18] and Lorentz-violating gravity [10, 11] 
by introducing a formalism to calculate the sensitivi¬ 
ties of stars (including pulsars) in generic theories. As 
an application, we calculate them in the most general 
scalar-tensor theories (“Horndeski theories” or “gener¬ 
alized galileons”) [19-21] that have second-order field 
equations and are invariant under a shift of the scalar 
field (j), i.e. (f> ^ (j> -\- constant. These theories have re¬ 
ceived much attention because some of them provide a 
screening mechanism (“Vainshtein mechanism” [22, 23]; 
see also Refs. [24, 25]), which may permit modifying grav¬ 
ity on cosmological scales (possibly reproducing cosmo¬ 
logical data without Dark Energy) while recovering GR 
on small scales, where gravitational modifications would 
be “screened”. Also, galileon interactions arise in the 
decoupling limit of massive gravity [26, 27]. 

We show that in shift-symmetric Horndeski theories 
(SSHTs) stellar sensitivities vanish^, and the leading- 
order GW emission matches GR’s. We conclude that 
SSHTs reproduce existing binary-pulsar data (provided 
that a PN expansion over Minkowski space is adequate 
- i.e. the binary’s dynamics is perturbative -, and that 
cosmological-expansion effects can be neglected [29-32]), 
but deviations might still appear for the sources targeted 
by upcoming GW interferometers. 

We use Latin (Greek) lower-case letters for space 
(spacetime) indices, and capital Latin letters for indices 
running on fields. Repeated indices denote summations, 
and we assume c = 1 and signature (—, -I-, -I-, -I-). 

A general expression for the sensitivities: Gonsider an 
action 


S = J C{QA,d^QA)d^x, (2) 

where Qa are fields, and the corresponding field equa¬ 
tions 

A solution’s canonical mass-energy is 

M = y d^x {iTAdtQA - C) , (4) 


where -ka = d£/d{dtQA) for brevity. For a solution with 
TTAdtQA = 0 (this is the case e.g. for stationary solutions, 
where dtQA = 0, or solutions where some variables Qa 
depend on time, but tta = 0, see e.g. Eq. (8)), the mass 


is then simply^ 

M = -J d^xC{QA,df,QA). (5) 

Consider a neighboring solution, with TTAdtQA = 0 and 
the same total baryon number and entropy [c.f. Eq. (1)]. 
The two solutions represent the same star with different 
local values of the fields Qa, and their mass difference is 

5M = - j d^xdt (ttaSQa) - j d^S,-^^^5QA , ( 6 ) 

where 5Qa is the difference between the fields, dfSi is 
a coordinate surface element, and we have used Eq. (3) 
to show that the bulk terms vanish, as well as Gauss’s 
theorem. 

Because this result assumes an action with no deriva¬ 
tives higher than first order, it would not seem to apply 
to GR, since the Einstein-Hilbert action depends on sec¬ 
ond metric derivatives. However, the latter enter the 
Einstein-Hilbert action only through a total divergence, 
i.e. GR can be described by a first-order “Einstein La- 
grangian” Cg = -r^^r^„)/(167rG) (see 

e.g. Ref. [34]). Besides the metric, in theories different 
from GR there are other gravitational degrees of freedom, 
which we denote by (pA, and which we assume coupled 
to the metric and its derivatives, i.e. with Lagrangian 
C^{(l)A,dg,(t>A,gtiu,dagtj.,y). Moreover, the matter fields 
■05 must couple minimally to the metric [i.e. with La¬ 
grangian /lm(0B)ff/ijv)] to satisfy the weak EP^. 
Consider a stationary (i.e. time independent) star in one 
such theory. The mass difference between neighboring 
solutions is 


SM = - d^S, 


dCa 


d{d,g^,,) 


dgfj 


- J d^xdt (tt^bJ^b) - J d^St 


- / 




d{di^A) 


d{d^ipB) 

S^A- f d^S 5g^,, (7) 


where we have used dtijTg^^dggv) = dt{TT^^StpA) = 0 be¬ 
cause of stationarity. Note that we have not assumed 
dtiTTiPj^Stps) = 0 [but only TT^j^dttfB = 0, so that Eqs. (5) 
and (6) hold], which will allow using this expression for 
perfect-fluid stars. 

The Lagrangian for a perfect fluid with equation of 
state p = p(n, a) (p, n and a being respectively the 
energy density, baryon-number density and entropy per 
particle) is [35, 36] 


£771 = -0dg{aJ^)-aAdg,{/3AJ^), (8) 


^ This agrees with previous results for a specific theory of this 
class, i.e. Einstein-dilaton Gauss-Bonnet gravity in the decou¬ 
pling limit [28]. 


^ We have checked that this mass matches the ADM mass, both 
in GR and SSHTs (see also Ref. [33]). 

^ The fields (pA and 'ips are not necessarily scalars, e.g. they could 
represent the components of a vector or tensor. 
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where Pa, and B are scalars {a a can be interpreted 
as Lagrangian coordinates). The fluid four-velocity is 
defined as /\J\ (with |J| = yjJ'') and 

n = \ J\/^—g, i.e. is the baryon number density cur¬ 
rent Variation with respect to (keep¬ 

ing (fi, 6, a, a A, Pa fixed) gives the perfect-fluid stress- 
energy tensor, by using the first law of thermodynamics 
ndp/dn\a = p + p, {p being the pressure). Variations 
with respect to J^, ip, 9, a, aA, Pa yield 


= d^{cjJ^) = Jf^d^PA = J^d^aA = 0 , 
hU^ = -d^p - ad^e - P^d^aA , U^^^^B = T, (9) 

where h = {p + p)/n\s the specific enthalpy, and we have 
used dpjda\n = nT (T being the temperature) from the 
first law of thermodynamics. From these equations, it 
is clear that p and 6 are Lagrange multipliers enforc¬ 
ing the local baryon-number and entropy conservation, 
while Pa and a a are constant along the fluid lines. Also, 
these equations imply the conservation of the fluid stress- 
energy tensor [35, 36]. 

For a stationary fluid, we can adopt comoving coor¬ 
dinates where C/® = J® = 0. Therefore, the third term 
in Eq. (7) vanishes (since dJ® = 0), while the second 
becomes 


= J d^xdt[p6J*+ e5{aJ*) + aAS{P^J*)] 

= - j d^x [{h - aT)UtSJ* + TUtS{aJ*)] , (10) 

where we have used 9* a a = dtP^ = dtJ* = dtcr = 0, 
—dtp = {h — aT)Ut and —dtB = TUt, obtained from 
Eq. (9) (with J® = 0). Eor a fluid in hydrostatic and 
thermodynamic equilibrium in a stationary spacetime, 
TUt and {h — aT)Ut are uniform [37-39], and Eq. (10) 
thus becomes 

- J d^xdt = -ih-aT)UtSN-TUtST , ( 11 ) 


where 6N = J dPx5J^ and SS = f d^x6{aJ*) are the 
differences in total baryon number and total entropy be¬ 
tween the two solutions. Therefore, the terms in Eq. (7) 
depending on the matter variables vanish, if the two so¬ 
lutions have the same entropy and baryon number [c.f. 
Eq. (1)] 

Also, if the two solutions are asymptotically flat, 
i.e. g^iu = Ptiu + 0[\/r), the first term in Eq. (7) also 
vanishes [17]. This is seen by evaluating the integral 
at r —>• oo, since the integrand decays as 1 /r^, while 
d^Sr ^ Therefore, the mass variation only depends 
on the non-GR part of the action, i.e. 
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didipA) 



d^S, 


dC^ 


( 12 ) 


This generalizes similar expressions for scalar-tensor the¬ 
ories [17], and Lorentz-violating gravity [10, 11]. 

Let us now consider an action depending also on second 
derivatives of the metric and extra gravitational degrees 
of freedom pA (e.g. Horndeski theories): 

S = J C{QA,df,QA,di,df,QA)d'^x. (13) 

By introducing new fields = d^QA and enforcing 

this definition by Lagrangian multipliers, one obtains 

S= [ [C{Qa 5.A^^) + X^^iXAt. - d^QA)] dA , 

(14) 

whose variation relative to Xa^i yields 

Since Eq. (14) is in the form given by Eq. (2), the con¬ 
struction outlined above gives [using also Eq. (15)] 
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dPnv ■ (16) 


Sensitivities and gravitational radiation in SSHTs: 
SSHTs are described by the Lagrangian for the “galileon” 
scalar p [19-21]: 


Cd> = 


VEl 

IQttG 
-\-G 4 x 


[k{X) - Gii{X)Up + G4{X)R 


{upf - [X^X.pf 
G5X 


+G 5 (A)G^,V^V"(^- 


(□<^)^ 


-3 (□(/)) {Vt^V.pf + 2 (V^V,<)))"] + xPQ), (17) 


where V, R and G^i/ are the Levi-Civita connection, 
Ricci scalar and Einstein tensor, K, G 3 , G 4 , and G 5 
are arbitrary functions of A = — V^(/)V'® 0 / 2 , Gix = 
dG,/dX, □ = V^V^, (V^V,0)^ = 

(V^Vjyi)))^ = V^jS/PpVpV''pVv^^p, X is a constant, 
and Q = R^^Xk. — d:R^'"Rpu + R? is the Gauss- 

Bonnet scalar^. The total action also includes the Ein¬ 
stein Lagrangian Lg and the perfect-fluid Lagrangian 


The xdS term is shift-invariant because Q is (locally) a total 
divergence. Also, this term can be obtained by choosing G 5 oc 
ln|X| [40]. 


5X1 = - 
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Cm described above. Regarding the latter, one may 
couple the matter fields to a “disformal” metric g^^i, = 
being a constant), rather than to 
alone. If such a disformal coupling is present, however, 
one can adopt g^^, as the metric field, which puts the 
matter Lagrangian in the form Cm{'4’B,d^'4>B,g^iv) con¬ 
sidered above, while the action (17) remains invariant up 
to redefinitions of the functions K, G 3 , G 4 , and G 5 [41]. 
Thus, our results also apply to SSHTs with a (special) 
disformal coupling to matter. To allow asymptotically 
fiat solutions (see below), we assume that G 3 , G 4 , and 
G 5 are analytic in X This implies K{X) = X + 0{XY^ 
since a constant can be absorbed in the matter stress- 
energy tensor as an effective cosmological constant, and 
a coefficient for the linear term can be absorbed by re¬ 
defining (j)] Gz = 0{X) and G 5 = 0{X), since a constant 
produces a total divergence®; and G 4 = 0(X), since a 
constant in G 4 can be absorbed in the metric Lagrangian 
Cg by redefining the “bare” Newton constant G. 

Consider now an isolated stationary star, i.e. (j) = 
G(l/r), g^i, = r]fj,^+0{l/r), where the shift symmetry al¬ 
lows setting (j) asymptotically to zero. [In theories with a 
Vainshtein screening, cj) = 0{l/r) and = 77 ^j/-|-G(l/r) 
only for r ^ with the Vainshtein radius within 
which deviations from GR are screened; we will return 
to this later.] To determine the sensitivities, recall that 
Eq. (16) compares two neighboring solutions. If the lat¬ 
ter are asymptotically flat, ~ i? ^ ~ R^vap = 

G(l/r)®, and the differences 5(j) and Sg^y between them 
scale as 5(j) = 0(l/r) and Sg^y = G(l/r). With these 
asymptotics, all the surface integrals in Eq. (16) vanish 
when evaluated at r —> 00 , thus yielding zero sensitivi¬ 
ties. For example, by considering the contribution of K 
to the first term in Eq. (16), we find 

and this surface integral vanishes when evaluated at r —>• 
00 . Similar calculations show that all terms in Eq. (16) 
vanish. 

Note that it is the shift symmetry that makes the sen¬ 
sitivities vanish. In generic Horndeski theories, two solu¬ 
tions (“ 1 ” and “ 2 ”) have in general different asymptotic 
values of (j), i.e. = 4>ao + ot/r 0 (l/r^) (a being a 
constant) and = (j)oo + ^<^00 + («-!- 6a)/r + G(l/r^). 
As such, S(j) = S(l>oc + C){\/r) (while 5(() = 0{\/r) in 
SSHTs) and Eq. (18) gives 5M oc a5(j)oo, hence oc 
dM/d(j)oo oc a. This is e.g. the case for FJBD theory [5-7] 
and Damour-Esposito-Farfee gravity [18], where the sen- 


® This excludes e.g. K{X) ~ which reproduces Modified 

Newtonian Dynamics (MOND) [42] in the non-relativistic limit, 
and gives (j) \nr near spatial infinity. 

® Recall the Bianchi identity = 0. 


sitivities are proportional to the coefficient a of the l/r- 
term in the scalar field’s fall-off. Therefore, these theories 
predict the emission of dipolar radiation, hence they can 
be constrained by existing binary-pulsar data [12], and 
give testable predictions for upcoming GW interferome¬ 
ters [14-16, 43, 44]. Also, in general the sensitivities are 
not zero in the presence of a conformal coupling between 
the galileon and matter, e.g. in massive gravity. This 
effect is not considered e.g. in Refs. [45, 46]. 

The vanishing sensitivities imply the absence of dipo¬ 
lar gravitational emission in SSHTs. Gonsider metric and 
scalar perturbations hg_y and 5(j) over a Minkowski back¬ 
ground, i.e. g^y = rif,y+ehgy+0{e)^ and cj) = e6(j)+0{e)‘^, 
with e a perturbative parameter. To leading order in e, 
Eq. (17) coincides with the Lagrangian of a minimally 
coupled scalar field^, i.e. the leading-order field equations 
for a binary become 

Dnh^‘' = -16 ttGT>^^ + , (19) 

ngSct> = , sf) + eO{U(j>, 54>\ h^) 

= eO{h5cl>,6(j)^,h'^), ( 20 ) 

where = r]^'^dg,dy, — 

= 0 are the sensitivities, 

is the binary’s stress energy tensor (which is the same 
as in GR and depends on the stars’ masses and veloci¬ 
ties), and we have assumed the Lorenz gauge = 0. 

(These equations can also be obtained from Ref. [17], by 
noting that the theories studied here coincide at lead¬ 
ing order with those of Ref. [17], when the conformal 
scalar-matter coupling is switched off, i.e. A{(j)) = 1 in 
Ref. [17] ’s notation.) Equation (20) implies that (j) is not 
excited at leading order in e, while Eq. (19) matches its 
GR counterpart. As such, gravitational emission behaves 
as in GR at leading PN order, i.e. monopolar and dipo¬ 
lar emission vanish - like in GR, but unlike non-shift- 
symmetric scalar-tensor theories, where these emission 
channels are sourced by the sensitivities [17, 47] -, while 
the quadrupolar emission matches GR’s®. 

Replace now the leading-order solution in the non¬ 
linear terms of Eqs. (19) - (20). Since = 0 at leading 
order, at next-to-leading order one has 

= -167rG[(l - + er^""] -f G(e)2 , (21) 

= -exSg + Oief . ( 22 ) 

Here, = 0(h?) is GR’s gravitational stress-energy 
pseudotensor [48], while SQ = 0(h)^ is the perturbed 


^ Note that up to boundary terms, G 4 = 0{X) gives 0(e)^- 
terms in the action, because XR + = 

+total divergence [40]. 

® FJBD theory corrects GR’s quadrupole formula even for = 
= 0, due to the (Einstein-frame) conformal scalar-matter 
coupling. This coupling (and corresponding corrections) are not 
present here. 
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Gauss-Bonnet invariant. Equations (21) - (22) match 
those of Einstein-dilaton Gauss-Bonnet gravity [i.e. K = 
X and G 3 = G 4 = G 5 = 0 in Eq. (17)] in the decoupling 
limit, hence the leading PN order at which deviations 
from GR appear is the same as in that theory. Following 
Ref. [49], we then conclude that non-GR effects only ap¬ 
pear at 3PN (2PN) order in the dissipative (conservative) 
sector. 

Our PN formalism is only valid for stellar radii R 
much smaller than the gravitational wavelength , so 
that the perturbations decay as l/r in both the far zone 
(i.e. at distances from the binary r ^ A^^.^^) and near zone 
(i.e. for i? <C r A^,^) [50]. In theories with a Vain¬ 
shtein mechanism, (j) = 0(l/r) and = 77 ^ 1 / + 0(1/r) 
only for r ^ rv, i.e. is an “effective” stellar ra¬ 
dius. Therefore, our perturbative PN approach requires 
Tv <C in which case only causes higher-PN or¬ 

der “finite-size” effects [50]. Note that for the dominant 
quadrupole mode, A(,.^^ ^ 10 ® km for binary pulsars and 
Aqw ^ 10 ^ km for the late inspiral of neutron-star bina¬ 
ries targeted by upcoming GW detectors. Thus, although 
the value of (and the very presence of a Vainshtein 
mechanism) depend on the theory, our approach might 
break down, especially in the latter case. If the dynam¬ 
ics is not perturbative/PN the analysis may be conducted 
case by case using a WKB approach [45, 46]. While sim¬ 
ple theories (e.g. cubic galileons) may still provide results 
in agreement with binary-pulsar data [45] [see however 
also Ref. [51]], in more generic theories (e.g. ones includ¬ 
ing quartic and quintic galileons), the non-perturbative 
dynamics generally makes the WKB approach also fail 
unless rv A^^^v [46]. If this condition is not satis¬ 
fied, Ref. [46] concludes that many multipoles radiate 
with the same strength. This seems difficult to recon¬ 
cile with binary-pulsar observations, which agree with 
GR’s quadrupole formula. (Note however that unlike us, 
Ref. [46] assumes a conformal scalar-matter coupling.) 

In conclusion, GW emission from stellar binaries is 
only modified at high PN orders in SSHTs. There¬ 
fore, binary-pulsar observations, which agree with 
GR’s quadrupole formula at percent-level and also test 
the IPN conservative dynamics, are reproduced in 
SSHTs, provided that the binary’s dynamics is per¬ 
turbative/PN [46], and that cosmological-expansion ef¬ 
fects can be neglected [29-32]. Direct detection of GWs 
from neutron-star binaries, however, may still provide 
prospects for testing these theories. Also, our main re¬ 
sult Eq. (16) applies to generic gravitational theories, but 
only to stars and not black holes (where additional sur¬ 
face integrals at the horizon might be present). In the 
special case of SSHTs, black holes are the same as in 
GR if y = 0 in Eq. (17) [52], but possess scalar hairs if 
X 7 ^: 0 [53-55]. We thus expect dipolar emission from bi¬ 
naries involving black holes in SSHTs with x 0 [49, 56]. 
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